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Abstract. We introduce a new rule of motion for a totally asymmetric exclusion 
process (TASEP) representing pedestrian traffic on a lattice. Its characteristic feature 
is that the positions of the pedestrians, modeled as hard-core particles, are updated 
in a fixed predefined order, determined by a phase attached to each of them. We 
investigate this model analytically and by Monte Carlo simulation on a one-dimensional 
lattice with periodic boundary conditions. At a critical value of the particle density a 
transition occurs from a phase with 'free flow' to one with 'jammed flow'. We are able 
to analytically predict the current-density diagram for the infinite system and to find 
the scaling function that describes the finite size rounding at the transition point. 
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1. Introduction 

Let a set of hard-core particles, labeled by an index i = 1,2, . . . , N, move 
unidirectionally from site to site on a one-dimensional lattice. We imagine 
that the particles represent pedestrians all walking at the same pace but not 
necessarily in phase with each other. This leads us to formulate the following rule of 
motion, that we state as the update scheme of a Monte Carlo simulation. Each particle 
i is assigned a phase G [0, 1), permanently attached to it, and during each time step 
(that is, each unit time interval) all N particles make a forward hopping attempt in 
the order of increasing phases. An attempted hop will succeed only if the target site is 
empty. 

This model is an instance of what are commonly called Totally Asymmetric Simple 
Exclusion Processes (TASEP); its novelty resides in its update rule. Before continuing 
the discussion, we mention some connections to existing work. 

Processes of particles moving stochastically on - often one- dimensional - lattices 
serve on the one hand as archetypes of out-of-equilibrium systems, and on the other 
hand as modeling tools to study transport in various systems, ranging from road 
and pedestrian traffic to intracellular traffic [1]. The particle motion may take place 
according to a large diversity of hopping rules. By the 'exclusion' principle one imposes 
the hard core condition (at most one particle per site); the 'total asymmetry' forbids 
backward hops; and the process is called 'simple' when hops are only between nearest- 
neighbor sites. 

Given these three properties that are characteristic of a TASEP, it is still possible 
to choose from a variety of update schemes. In particular, the following update schemes 
have been studied: parallel update [21 El IH E], random sequential update, sequential 
update ordered backward or forward in space [HI [71 El 13 13] , sublattice update [IDl El [I2] , 
and random shuffie update [131 HI]. The properties of the system depend on the update 
scheme [15] and the choice of the scheme should be determined by the application. 

The most common update schemes are the random sequential and the parallel 
updates. Random sequential update produces a dynamics very close to that defined 
by a master equation in continuous time. A time step is defined as a succession of 
N elementary updates, each associated with a time interval of length 1/A^, and each 
allowing only a single particle, chosen at random, to make a hopping attempt. With 
this dynamics considerable fiuctuations occur, since the same particle may be updated 
several times in the same time step whereas another one may be ignored during several 
time steps. 

With parallel update particles make hopping attempts only at integer values of time 
but then do so simultaneously. Parallel update is used in particular for applications to 
road traffic [T6l [T7] : all vehicles are moving at the same time and the time step of 
the scheme is then supposed to represent a reaction time. Fluctuations are reduced, 
but parallel update can create conflicts - that should be settled by additional rules - 
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when more than one particle tries to hop onto the same target site. This may occur 
in particular in applications to pedestrian traffic, which usually takes place in two- 
dimensional space. 

In order to overcome the limitations of these two types of updates, the so-called 
'shuffie update' has been proposed for modeling pedestrian flow. In the random shuffle 
update [131 before each time step the particles are pre-arranged in a randomly 
chosen order and then each of them is updated once, exactly in that order. This update 
scheme was used for example in [18] for large-scale simulations of pedestrians. 

In the present paper we explore a variant of the shuffle update for which the order 
in which the particles are arranged, that is, the updating order, is fixed once and for 
allil. Our scheme is therefore appropriately characterized by the name of ^frozen shuffle 
update'. 

For a closed system this is easy to implement; a random phase Tj is drawn for 
each of the particles independently, for example from the uniform distribution on 
[0, 1). In each time step all particle positions are updated once, one after the other, 
according to increasing values of their phases. The phases Tj do not change during the 
whole simulation and may be considered as frozen variables of the motion. The set 
{xj} determines a random permutation of the particles; for uniformly distributed Ti , all 
permutations have the same probability. 

A closed system is expected to evolve towards a stationary state. We must be 
prepared to envisage that the final stationary state may depend (and as we shall see, 
indeed does depend) on the precise permutation that fixes the updating order of the 
particles. An average over all permutations is therefore appropriate and is analogous to 
the averages on quenched disorder variables standardly performed in statistical physics. 
The term 'disorder average' will therefore denote below the average over all random 
assignments {rj}. 

For an open system the frozen shuffle update requires that by a suitable algorithm 
we fix the phase of each particle the moment it enters. The equivalence of the set {ti} to 
a simple permutation may then no longer hold. The case of open boundary conditions 
will not be considered here but is studied in a forthcoming paper [22] . 

In section [2] we introduce some terminology that actually already is the expression 
of several model properties. In section [3] we consider the TASEP with frozen shuffle 
update on a ring with particle density p. We show that a phase transition occurs at 
a critical density pc which separates a low density regime with 'free flow' from a high 
density regime with 'jammed flow'. We determine the current- fersus-density curve Jl{p) 
analytically, fflst for an infinite system (section [311!) where J(p) = lim^^-^oo Jl{.p) has a 
cusp, and then for a system of finite size L (section 13.21) . where the finite size rounding 
of Jl(p) is described by a scaling function that depends only on the product variable 
(p — Pc)L^^^ ■ Monte Carlo simulations show very good agreement with theory. In section 

I This frozen variant was mentioned, but not studied, in the conclusion of |13j . 
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Figure 1. Lattice sites are represented by squares that may be either empty or 
occupied by a single particle. A configuration involving four particles is shown 
at two successive times t = s and t = s + l; the flow is in the direction of the 
arrow. Particles 1 and 2 form a well-ordered pair. They can move at the same 
time step, as particle 1 is updated before particle 2. Particles 3 and 4 form 
an ill-ordered pair. During the time step from t = siot = s + l the update 
attempt of particle 4 is performed at time s + ta but remains unsuccessful, since 
its target site is still occupied by particle 3. When subsequently particle 3 is 
updated at time s + T3, it moves forward. Thus a hole is inserted between the 
two particles. In case rs < T2, particles 2 and 3 also form an ill-ordered pair, 
but they are not adjacent and thus can hop independently. 



m by way of a supplement, we show that under 'free flow' conditions the TASEP with 
frozen shuffle update is equivalent to a system of noninteracting particles in continuous 
space and time. In section |5] we conclude. 

2. Free flow and jammed flow 

We introduce here the concepts that characterize the structures formed by the particles 
as they result from the frozen shuffle update scheme. The discussion below is 
independent of any boundary conditions that may be imposed at a later stage. The 
most important point is the identification of two distinct stable flow states that we call 
the free flow state and the jammed state. 

2.1. Well-ordered and ill- ordered pairs 

Let the flow direction be to the right, let the particles be numbered . . . , — 1, i + 1, . . . 
from right to left (see flgure [1]), and let their phases . . . , rj_i, Tj, Xj+i, ... be given. The 
pair of successive particles (^,2 + 1), not necessarily on adjacent sites, will be called well- 
ordered if Tj+i > Ti and ill-ordered in the opposite case. The time evolution of well- and 
ill-ordered pairs under the frozen shuffle update scheme has the following properties, 
illustrated in figure [H 

If a well-ordered pair + 1) occupies two adjacent sites, then at each time step 
particle i will move first and particle i + 1 will move next; hence when the time step is 
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completed, the two particles are still adjacent and have advanced one lattice distance 
to the right. Their speed is f = 1 in units of lattice distances per time step. 

If an ill-ordered pair occupies two adjacent sites, the two particles cannot move in 
the same time step; particle i + 1, having Tj+i < Xj, will attempt first to move but finds 
itself blocked by particle i. Hence the two particles of an ill-ordered pair move at speed 
V = 1 only if they are separated by at least one empty site. 

2.2. Free flow configuration 

A particle configuration will be said to satisfy the free flow (FF) condition when each 
ill-ordered pair has its two members separated by at least one hole. In view of the above, 
such a configuration is again identical to itself at the end of each time step except for 
a translation by one lattice distance to the right. This corresponds to a free flow with 
speed f = 1, and hence for a FF configuration we have 

J = P, (1) 

where J is the current and p the particle density. It is tacitly understood here that 
these quantities refer to time averages in a stationary state. 

2.3. Rising sequences and platoons 

Going along the lattice from right to left one may divide the particles encountered 
into sequences of increasing phases (for short: rising sequences). The set of particles 
(i, i + 1, . . . , i') will be said to constitute a rising sequence if < Tj+i < . . . < r,' but 
rj_i > Tj and Tj/ > rj/_|_i. Examples are shown in figure [2l 

Let a set of particles occupy consecutive sites and have phases that increase from 
right to left. If this set corresponds to a full rising sequence, it will be called a platoon. 
If it corresponds to only part of a rising sequence, it will be called a suhplatoon. One 
may say that a platoon (a subplatoon) is a fully compacted rising sequence (part of a 
rising sequence). Platoons and subplatoons are limited on both ends either by holes 
or by ill-ordered pairs. A rising sequence is composed either of a single platoon or of 
several subplatoons, an isolated particle being considered as a (sub-)platoon of length 
1. 

Under the frozen shuffle update platoons and subplatoons obey the following simple 
rules. 

(i) If in a given time step the first particle of a (sub-)platoon can move, then all 
its other particles will also move; hence (sub-)platoons move as single entities. 

(ii) When two subplatoons merge, they can never separate again; hence the length 
of a subplatoon can only grow until it includes the whole rising sequence in which it is 
embedded. 
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Figure 2. A six particle configuration is shown at two successive times t = s and 
t = s + 1. Because of tlie inequalities between their phases, the set of particles 
(1,2,3) forms a rising sequence, and so does (4,5,6). The inequality T3 > 
defines the separation between these two sequences. The fact that particles 3 
and 6 are the last ones of their rising sequences is marked by a heavy (red) 
line segment delimiting their lattice site to the left. At the {s + l)th time step 
all particles will move except number 4. It so happens that at time t = s + 1 
particles 4, 5, and 6 have formed a platoon, i.e. the rising sequence (4, 5, 6) has 
been compacted. 
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Figure 3. A jammed configuration involving seven particles is shown at two 
successive time steps t = s and t = s + 1. All particles are grouped together 
in platoons; particles belonging to the same platoon are labeled by the same 
letter. The last particle of each platoon is indicated by a heavy (red) line 
segment to its left. Successive platoons are separated by either zero or a single 
hole, that is, the particles are in a jammed configuration. During the (s + l)th 
time step platoons a and b move one lattice distance to the right, but c and 
d are blocked. Inversely, one may describe this dynamics as a motion of holes 
that jump at each time step across the platoon to their left. 
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2.4- Jammed configuration 

A configuration of particles will be called jammed if all its rising sequences are compacted 
into platoons and if consecutive platoons are separated by at most one hole. Figure [3] 
shows an example of a jammed configuration. At each time step the evolution of a 
jammed configuration may be simply described in terms of the motion of its platoons. 
The rules follow directly from those above: 

(i) A platoon preceded by a hole advances by one site as a single entity; this 
amounts to a position exchange of the platoon and the hole. 

(ii) A platoon not preceded by a hole is blocked and does not advance. 

Let V stand for the average platoon length in a jammed configuration that is 
statistically homogeneous in space. Noting that 1 — p is the hole density and that only 
platoons preceded by a hole move, we may write the particle current in this jammed 
configuration as 



Obviously, configuration space has many configurations that are neither 'free fiow' 
nor 'jammed' in the sense of the above definitions. There is also one subclass of 
configurations that are both 'free fiow' and 'jammed'; this happens when all platoons 
of a jammed configuration are separated by exactly one hole. 

3. Phase transition on a ring 

After the preliminaries of section |2] we are now ready to study a concrete system. We 
consider a ring, that is, a lattice of L sites with periodic boundary conditions. Let be 
the number of particles and hence p = N/ L the particle density. We set ourselves the 
purpose of determining the particle current Jl{p) as a function of the particle density 
p in the stationary state that will result from a given initial state. 

At the initial time t = the particles are placed at distinct but otherwise random 
positions on the lattice. They are numbered i = 1,2, . . . , N from right to left (clockwise 
around the lattice) and their direction of motion is from left to right (anticlockwise). The 
particles are assigned a random phases which we take independently and uniformly 
distributed on [0, 1). This assignment determines the updating order of the particles. 

The initial configuration does not necessarily satisfy the FF condition. If it does, 
then the particle configuration at time t = s is obtained from that at time t = by 
rotating all particle positions by s steps along the ring. If the FF condition is not 
satisfied, then after a transient period the system will reach a stationary state which 
may or may not be of the FF type. We will investigate below the conditions for the 
realization of each of these possibilities, and the ensuing consequences for the particle 
current. 




(2) 
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3.1. Infinite system limit 

The infinite system limit is easiest to discuss, since we may apply the law of large 
numbers and formulate statements that in that limit are true with probability 1. Let 
us first ask up to which value of the density p it is still possible to have free flow. 

For a given set {tj} the densest possible FF configuration occurs when all rising 
sequences are compacted into platoons separated by a single hole. This corresponds 
precisely to the aforementioned special case of a configuration which is both FF and 
jammed. The maximum density of the FF phase thus is pc = z^/ (z/+l), where as before i' 



is the average platoon length. It may be shown (see [19] or Appendix Appendix A ) that 
in the infinite system limit one has i/ = 2 when the phases Tj are uniformly distributed, 
and therefore pc = | • 

For p < Pc any arbitrary initial configuration - tacitly understood to be statistically 
homogeneous in space - will, after a transient, be converted into a FF configuration. 
Indeed, whenever an ill-ordered pair of particles occupies two successive sites, the second 
one will not yet be able to move when the first one first moves forward, and a hole will 
naturally be included between them. When in this way all ill-ordered pairs have come 
to include a hole, a FF configuration is obtained. The current J(p) = limi^oo Jl{p) is 
then given by its FF value ([1]), 

J{p)=P, P<Pc- (3) 

For p > Pc the time evolution will produce two effects. It will compact rising 
sequences into platoons and it will distribute the available holes such that each platoon 
is separated from its predecessor by at most a single hole. However, the number of holes 
is less than the number of platoons. The number of platoons that move in a given time 
step has thus been maximized and is equal to the number of holes, the other platoons 
being blocked at that time step. This corresponds to the definition of a jammed phase 
given in section 12. 4[ whence upon applying ([2]) with u = 2 we obtain 

J(p) = 2(1 -p), p>pc. (4) 

Equations ([3]) and (jll) lead to the cusped current-density diagram shown in figure HI 
The agreement with finite size Monte Carlo simulations is already quite good for system 
size L = 12. However, finite size effects are visible around the maximum, as shown in 
the inset of figure HI In the next section we shall refine the theory to account for this 
rounding of the transition. 



3.2. Finite system 

We consider in this subsection a finite ring of size L containing exactly particles; 
throughout we set p = N/L. Our interest is in the density dependent particle current, 
which in this finite system we shall denote by Jl{p) and whose definition we shall render 
precise. 

By the mechanism described above the system, whatever its initial configuration, 
will evolve so as to maximize the number of ill-ordered pairs that include a hole. For 
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Figure 4. Current J as a function of the density p for periodic boundary 
conditions. Solid line: theoretical prediction for an infinite system. Data points: 
Monte Carlo simulations for systems of size L = 12, 52, and 102. The inset is 
a zoom around the maximum of the curve. 



densities p < 1/2, there is enough space in the system to place a hole between each pair 
of particles. Then the FF condition can be fulfilled with certainty and the stationary 
state is a FF state. Denoting the current in the stationary state by Jnl we have 

^1 

JnL = J^= P, P<2- (5) 

For densities p > 1/2 it may or may not be possible to converge towards an FF 
configuration, depending on the random assignment {rj}. The considerations of section 
12.21 show that the discriminating quantity is the number of ill-ordered particle pairs in 
the initial state. We denote by ra" (by n^) the number of well-ordered (ill-ordered) pairs, 
so that 71^ + 11^ = N. It will be convenient to work with the difference variable 

n{{n})=n-~n\ (6) 

of which we shall henceforth suppress the argument. Because of the periodic boundary 
conditions, there is always at least one ill-ordered pair and one well-ordered pair in the 
system, so that n may take the values = 1, 2, . . . , — 1. A necessary and sufficient 
condition to fulfill the FF condition in the stationary state is to have at least one empty 
site available for each ill-ordered pair, that is, -|- 2n' < L or equivalently 

n>lN-L. (7) 

The expression for the stationary state current now involves the variable n and we will 
denote it by Jnlu- Two cases have to be distinguished. First, if inequality ([7]) is satisfied. 
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the system evolves towards a FF state and for this subset of reahzations the current is 

Secondly, we consider realizations {tj} for which inequality ([7]) is violated. The 
stationary state then only has isolated holefl and all rising sequences are compacted 
into platoons. At each time step only the platoons headed by one of the L — N holes 
move forward, which means that the instantaneous current per time step fluctuates with 
time. However, averaged over time each platoon will move in a fraction {L — N)/N of all 
time steps. Using the fact that N/n^ = 2/(1 — n/N) is the average length of a platoon, 
we therefore find after time averaging for the current Jnlu the expression 

L~N N 2(1- N/L) ^ n ^ -, 

L 1 — n/N ^ i\ ^ 

For finite systems and for densities | < p < 1, there will always exist realizations of the 
Ti that converge towards FF stationary states with a current p, and others that do not 
and have a current less than p and given by ([9]). In this density regime we will denote 
by Jnl the current Jnlti of ([8]) and ([9]) averaged with respect to n, that is, 

N-l 

Jnl = PN{n)JNLn , P>i> (10) 

n=l 

in which P/v(n) is the probability distribution of ^({tj}) and remains to be determined. 
Since n is determined by {rj}, the current J^l in ( !T0|) deserves the name of 'disorder 
averaged current'. 



3.3. Finite size effects near the transition point 

The probability distribution PN{n) was studied by Oshanin and Voituriez [19] for the 
case - which is also ours - where the are drawn independently from a uniform 
distribution on [0,1). These authors showed, among other things, that in the limit 
of large and with n scaling as ~ N^^'^ the variable x = n/N^/"^ has the probability 
distribution 

n(x) = (3/27r)i/2exp(-|x2). (11) 
It is symmetric in n, as dictated by the left-right symmetry of the phase assignment. In 



Appendix [Appendix A| we derive equation (JTTil in a more direct way. 

From here on we shall consider the equations of the preceding subsection in the 
limit of large but finite A^, n and L, and fixed ratios p = N/L and x = n/N^/'^. We will 
conform to usage and take the system size L, rather than A^, as the independent large 
variable. In the limit in question we shall write Jnlu = Jl{p,x) and Jnl = Jl{p)- We 



may then reexpress the disorder averaged current (JTOl) as 

POO 

Jl{p) = / dxIi{x)JL{p,x). (12) 



§ The same behavior appears with random shuffle update. 
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IS, 



X I l-z(pL)-i/2' X<Xc{p), 

Jl{,P,x) = { > (13) 

p, X> Xc(p), 



in which 



a;,(p) = (pL)V2(p-i_p-i) (14) 

where pc = |. We observe parenthetically that in the limit L ^ oo the x dependence 
of f|T^ disappears and we recover lim^^oo <^L(a;, p) = J{p), where J(p) is the infinite 
system current of equations and (jl]). Since f lTT]) is valid in the limit in which x 
remains finite as L — )■ oo, we conclude that the present approach is valid for densities 

p = Pc + Ap (15) 

such that Ap is on the scale of Remembering this and expanding in powers of 

we find from ([HD and ([13]) 

xcip) =(|)^/^LV^Ap + 0(L-V2), 

f pc-2Ap + (3L/2)-^/^x + 0(L-^), x < xJp), 
Mp,x)= , (16) 

{ Pc + Ap, X > Xc(p). 

We introduce the scaling variable y = L^^^Ap, which in the limit of interest should be 
of order unity. Substitution of f lT6|) in f lT2|) then yields 

JL{p)=Pc + L~'/'y-L-'/' / dxU{x) [3y-{iy/'x]+OiL-').il7) 

J — oo 

When using in (fTTl) the explicit expression ( ITT]) for n(x) we may evaluate the x integral 
and obtain, up to corrections of higher order in 

Jl{p)=P. + L-^"^{L^"Ap), (18) 

valid in the limits Ap = p — pc — ?■ and L — )■ oo with L^^'^Ap fixed, and in which the 
scaling function $(?/) is given by 

^{y) = -\y - fl/erf (fy) - (97r)-i/2exp (-f y^) . (ig) 
This function is negative and such that 

^{y) -y-, y -oo, 

<^{y) - -2y, y^oo, (20) 

which ensures the correct limit behavior of ( |T8]) for |p — pd ^ We have plotted 

$(?/) in figure [5] together with simulation data for different system sizes L. The data 
are seen to collapse very well onto the theoretical curve. 




Figure 5. Solid line: the theoretical scaling function ^{{p — pc)L^^^) of equation 
, representing the average current Jl {p) — Jc as a function of the particle 
density p in a finite system of size L near criticality. The dashed lines are the 
asymptotes for {p— pc)L^^'^ ±00. Simulation data for large system sizes L are 
seen to collapse very well onto the theoretical curve. Each point corresponds 
to an average over 10000 or 100 000 time steps and over 1000 realizations of 
the disorder. 



4. Mapping to a continuous model and interpretation for pedestrian motion 

In this section we point out that under free flow conditions the time evolution defined by 
the frozen shuffle update for the particle system on a lattice may be seen as a sequence 
of snapshots taken at integer instants of time t = . . . , s — 1, s, s + 1, . . ., of a system that 
itself evolves in continuous time t and space x. 

To show this we consider a collection of nonoverlapping hard rods all moving 
continuously to the right at speed v = 1 along the depicted in figure Ei 

If we associate lattice sites with the integer axis positions a; = . . . , — 1, fc, A; + 1, . . ., 
then at any given instant of continuous time, each rod covers exactly one site. The 
mapping is performed by placing on that site a particle associated with that rod. Let 
figure [6] represent the rod positions at time t = (or for that matter at any other integer 
instant of time). The particle labeled i and corresponding to rod i occupies lattice site 
k and therefore gives rise, at t = 0, to a particle on site k. The distance between site k 
and the tail of rod i has been indicated as a time interval]]] Xj, this being the time still 
needed for the tail of that rod to cross the point k during its continuous motion along 
the X axis. This crossing therefore occurs at time t = Ti, and that is the time at which 

II Because w = 1, times and distances may be identified. 
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Figure 6. Hard rods labeled by an index i move at constant speed v = \ along 
the X axis. The integer axis positions have been labeled by an index k. 



particle i will hop from site k to site fc + 1. Particle i will execute its subsequent hops at 
times t = s + Tj, where s is an integer. This is exactly the frozen shuffle update scheme. 

We remark that the mapping defined here yields only the FF configurations of 
the discrete model. If we try to perform the inverse mapping, i.e. from the discrete 
to the continuous model, then in case of a jammed configuration the non-overlapping 
condition for rods cannot be enforced anymore. This may actually still have some 
physical relevance, if one adopts the view that a rod represents not only a pedestrian 
but also some "private" space around him. In free flow pedestrians are not willing to 
approach each other too closely and they avoid to enter each other's "private" space, 
whereas at increasing densities they tolerate smaller distances. 

5. Conclusion 

We have introduced in this paper a new update scheme for the TASEP, namely the 
frozen shuffle update, which should be appropriate, in particular, for the modeling of 
pedestrians. 

We have characterized the behavior of the TASEP with frozen shuffle update for 
a closed one- dimensional lattice of L sites and N particles. The time evolution under 
frozen shuffle update is deterministicj^l; it is fully determined by the initial particle 
positions and by the set {tj} of their phases. The latter are quenched random variables 
that at each time step determine the update order of the particles. We showed that 
the analysis of the particle motion and their interaction may be fruitfully carried out in 
terms of the concepts of well/ill-ordered pairs and of platoons. Two principal types of 
flow may then be distinguished, 'free flow' and 'jammed flow'. 

We were able to predict completely the fundamental diagram, that is, the current 
Jl(p) as a function of density p = N/L, both for the inflnite (A^, L — )■ oo) and the flnite 
system. 

We found that for increasing particle density p the passage from a free flow phase 
to a jammed phase takes place via a phase transition at a critical density p = p^. 
This contrasts with the random sequential update, for which with increasing density 



^ By this we mean that the hopping probabihty is always unity when the target site is empty. 
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the system becomes gradually more and more congested. Critical points, however, were 
observed in the fundamental diagram in other instances of deterministic motion, namely 
with parallel update jSHlEl] and with random shuffle update [131 E]. the latter case, 
although the particle-hole symmetry is broken, the critical point was still found at the 
symmetric point pc = 1/2; by contrast, for the present frozen shuffle update we find 
Pc = 2/3, i.e. the asymmetry between holes and particles is still enhanced. Another 
difference is that for the random shuffle update the critical point is already present 
in finite systems, whereas for the frozen shuffle update the transition is rounded and 
becomes sharp only in the limit of infinite system size. 

A mapping with a continuous model of hard rods is proposed, which is exact for 
free fiow configurations, and may be useful for the interpretation of the results in terms 
of pedestrian motion. 

Two final remarks about open questions are in place here. First, the deterministic 
time evolution studied in this paper entails that, if the target site is empty, particles hop 
with probability p = 1. Whereas in the case of a random sequential update the hopping 
probability p can be modified through a simple rescaling of time, here such a rescaling is 
not possible. We therefore expect a qualitatively different behavior of the current J(p) 
when the hopping probability p is strictly less than one. We leave the analysis of this 
case for future work. 

Second, this work has been exclusively concerned with a closed system. New 
types of questions arise when one applies frozen shuffle update to open systems. In 
a companion paper [22] we shall address the case of open boundary conditions and 
determine in particular the phase diagram. 

Appendix A. Random walk generated by a random permutation of N 



We arrange the integers 1, 2, 3, . . . , on the sites of a circular lattice and permute them 
randomly, all permutations having the same probability. Suppose that when going 
clockwise along the lattice in A^ steps, we encounter well-ordered and ill-ordered 
pairs in the sense of section 12. 1[ Obviously and n' are random integers that depend 
on the permutation, and are such that ra™ + = A^. Let n = n" — n'. We ask what the 
probability distribution PN{n) of n is in the limit of large A^. 

This question was first asked by Oshanin and Voituriez [19], who obtained the 
distribution n(x) given in ( JTTl) . It is possible to arrive at same result in a different and, 
we believe, simpler way that we present here. It is based on establishing a recursion 
in A^. Suppose that the integers 1,2, ...,A^ have been permuted and placed on the 
sites of a circular A^-site lattice. A permutation of 1, 2, . . . , A^ -|- 1 on an (A^ -|- l)-site 
lattice is obtained by inserting between two randomly chosen neighboring sites a new 
site carrying the integer A^ + 1. The probability pj^ (or p\r) to perform the insertion on 
a well-ordered (or on an ill-ordered) pair is 



integers 
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In either case the original pair disappears and, since the newly inserted integer + 1 is 
necessarily larger than its two neighbors, is replaced with the succession of a well- and 
an ill-ordered pair. Hence we have the recursion 

PN+i{n) = p^^{n - l)P^{n - 1) + p'^{n + l)P^{n + 1), (A.2) 

valid for n = —N + 1, — A^ + 3, . . . , A^ — 1 (which are the only values of n that can occur) 
and with the convention that Pn{-N) = Pn{N) = 0. We substitute (lA.ll) in (1A.2I) and 

set 

^ = ]7v^ ' = mr2^N j , (A.3) 

expecting that in the large- A^ limit the variables x and A^ may be treated as continuous. 
On the expression thus obtained we perform a standard expansion in negative powers 
of A^. The result is the Fokker-Planck equation 

dUNjx) _3 dxIlNix) I g^njv(x) 

dN ~ 2 dx 2 dx^ ' ^ ■ ^ 

of which flTT]) is the stationary solution, that is, the one solving dIlN{x)/dN = 0. 

We also note that the average length of the platoons 

N 2N , 
u= — = A.5 

n' N-n ^ ^ 

tends to V = 2 when A^ becomes large, as n typically scales as N^^"^. 



[1 

[2] 
[3] 
[4] 
[5] 
[6] 

[7] 



[9 

[lo: 

[11 
[12: 



D. Chowdhury, A. Schadschneider, and K. Nishinari. Physics of transport and traffic phenomena 
in biology: from molecular motors and cells to organisms. Physics of Life Reviews, 2:318-352, 
2005. 

H. Hinrichsen. Matrix product ground states for exclusion processes with parallel dynamics. J. 

Phys. A: Math. Gen., 29:3659-3667, 1996. 
M.R. Evans. Exact steady states of disordered hopping particle models with parallel and ordered 

sequential dynamics. J. Phys. A: Math. Gen., 30:5669-5685, 1997. 
A. Schadschneider and M. Schreckenberg. Garden of Eden states in traffic models. J. Phys. A: 

Math. Gen., 31:L225-L231, 1998. 
M. Wolki and M. Schreckenberg. Exact matrix-product states for parallel dynamics: open 

boundaries and excess mass on the ring. J. Stat. Mech., P05014, 2009. 
A. Borodin, P.L. Ferrari, and M. Prahofer. Fluctuations in the discrete TASEP with periodic initial 

configurations and the Airy process. International Mathematics Research Papers, rpm002, 2007. 
J.G. Brankov, V.B. Priezzhev, and R.V. Shelest. Generalized determinant solution of the discrete- 
time totally asymmetric exclusion process and zero-range process. Phys. Rev. E, 69:066136, 

2004. 

J.G. Brankov, V.V. Papoyan, V.S. Poghosyan, and V.B. Priezzhev. The totally asymmetric 
exclusion process on a ring: Exact relaxation dynamics and associated model of clustering 
transition. Physica A, 368:471-480, 2006. 

V.S. Poghosyan and V.B. Priezzhev. The relaxation dynamics of the TASEP with particle- 
dependent hopping probabilities on a ring. Reports on Mathematical Physics, 61:239-246, 2008. 

V. Fayaz, F.H. Jafarpour, S.R. Masharian, and S. Zeraati. One-transit paths and the steady state 
of a non-equilibrium process in a discrete-time update. J. Stat. Mech., P12009, 2010. 

C. Pigorsch and G.M. Schiitz. Shocks in the asymmetric simple exclusion process in a discrete-time 
update. J. Phys. A: Math. Gen., 33:7919-7933, 2000. 

S.S. Poghosyan, V.B. Priezzhev, and G.M. Schutz. Green functions for the TASEP with sublattice 
parallel update. J. Stat. Mech., P04022, 2010. 



Frozen shuffle update for TASEP 



16 



[13] M. Wolki, A. Schadschneider, and M. Schrcckcnbcrg. Asymmetric exclusion processes with shufHed 

dynamics. J. Phys. A: Math. Gen., 39:33-44, 2006. 
[14] D.A. Smith and R.E. Wilson. Dynamical pair approximation for cellular automata with shuffle 

update. J. Phys. A: Math. Theor., 40(ll):2651-2664, 2007. 
[15] N. Rajewsky, L. Santen, A. Schadschneider, and M. Schreckenberg. The asymmetric exclusion 

process: Comparison of update procedures. J. Stat. Phys., 92:151-194, 1998. 
[16] K. Nagel and M. Schreckenberg. A cellular automaton model for freeway traffic. J. Phys. I, 

2:2221-2229, 1992. 

[17] C. Appert and L. Santen. Boundary induced phase transitions in driven lattice gases with meta- 

stable states. Phys. Rev. Lett., 86:2498-2501, 2001. 
[18] H. Kliipfel. The simulation of crowds at very large events. In A. Schadschneider, T. Poschel, 

R. Kuhne, M. Schreckenberg, and D.E. Wolf, editors. Traffic and Granular Flow '05, pages 

341-346, 2007. 

[19] G. Oshanin and R. Voituriez. Random walk generated by random permutations of {1, 2, 3, . . . , n + 

1}. J. Phys. A: Math. Gen., 37:6221 6241, 2004. 
[20] A. Schadschneider and M. Schreckenberg. Cellular automaton models and traffic flow. J. Phys. 

A: Math. Gen., 26:L679-683, 1993. 
[21] B. Eisenblatter, L. Santen, A. Schadschneider, and M. Schreckenberg. Jamming transition in a 

cellular automaton model for traffic flow. Phys. Rev. E, 57:1309 1314, 1998. 
[22] C. Appert-Rolland, J. Cividini, and H.J. Hilhorst. Frozen shuffie update for an asymmetric 

exclusion process with open boundary conditions. Preprint, 2011. 



